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The relation between bulk topological invariants and experimentally observable physical quantities 
is a fundamental property of topological insulators and superconductors. In the case of chiral sym- 
metric systems in odd spatial dimensions such as time- reversal invariant topological superconductors 
and topological insulators with sublattice symmetry, this relation has not been well understood. We 
clarify that the winding number which characterizes the bulk Z non-triviality of these systems can 
appear in electromagnetic and thermal responses in a certain class of heterostructure systems. It 
is also found that the Z non-triviality can be detected in the bulk "chiral polarization", which is 
induced by magnetoelectric effects. 

PACS numbers: 



An important feature of topological insulators (TIs) 
and topological superconductors (TSCs) is that topolog- 
ical invariants characterizing the bulk states emerge as 
physical quantities probed by electromagnetic or thermal 
responses For instance, the Chern number appears 
as the quantized Hall conductivity in the quantum Hall 
effect state Q , and the Z 2 invariant of a time- reversal in- 
variant (TRI) TI in three dimensions can be detected in 
axion electromagnetic responses Q . The correspondence 
between bulk topological invariants and electromagnetic 
(or thermal) responses naturally arises from the existence 
of underlying low-energy effective topological field theo- 
ries [H . For most classes of TIs and TSCs @, [BJ , this cor- 
respondence has been well clarified so far. However, for 
the case of TIs and TSCs characterized by Z invariants in 
odd spatial dimensions, this point has not yet been fully 
understood. These classes include time-reversal symme- 
try broken (TRB) TIs with sublattice symmetry in one 
and three dimensions (class AIII), TRI TSCs in three di- 
mension (class Dili, e.g. 3 He, C Ua; Bi 2 Se3 0,0, Li 2 Pt 3 B 
[1]), and TRI TIs and TSC of spinless fermions in one 
dimension (class BDI, e.g. Su-Schrieffer-Heeger model 
Kitaev Majorana chain model [HJ). It is noted that 
all of these classes possess chiral symmetry (sublattice 
symmetry); i.e. the Hamiltonian H satisfied the relation 
THT = —H with T a unitary operator. This implies that 
if is an eigen state of H with an energy E, then, T\ip) 
is also an eigen state with an energy —E. The chiral 
symmetry is indeed the origin of the bulk Z topological 
invariant referred to as the winding number. The chiral 
symmetric topological insulator with the winding num- 
ber N possesses N flavors of gapless Dirac (Majorana) 
fermions at the boundary, which are stable against dis- 
order and interactions, as long as the chiral symmetry 
is preserved 11] . To this date, however, it has not been 



case of TRI Z 2 TIs, the action of which is given by [3, EH , 



fully elucidated how the winding number can be detected 
in electromagnetic or thermal responses. For instance, in 
the case of three-dimensional (3D) class AIII TIs, low- 
energy effective theory is the axion field theory as in the 
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is the magnetoelectric polarization expressed by the 
Chcrn-Simons 3-form with Berry connection A nm {k) = 
i (u n (k) I du m (k)) for occupied states \u n (k)). Because of 
chiral symmetry, P3 takes only two values, i.e. P3 = 4^ 
(mod 1) where N3 is the integer- valued winding num- 
ber [ljj]. Thus, the above field theory captures only Z 2 
part of the winding number, and fails to describe the Z 
nontrivial character [15| . The same problem also occurs 



for class Dili TSCs, as previously noticed by Wang and 
his coworkers [lij]. For this case, Wang et al. presented 
an argument based on an effective theory for surface Ma- 
jorana fermions. However, a general framework which 
relates the winding number to electromagnetic or ther- 
mal responses is still lacking, and desired. In this paper, 
we present two approaches for the solution of this issue. 
One is based on the idea that the winding number can be 
detected in electromagnetic and thermal responses of a 
certain class of heterostructure systems (see Eqs.([7])-(|9]) 
and (TlT|l - (ri3|) below, which constitute main new results). 
We clarify the condition for the heterostructure systems 
in which the Z non-trivial character of the bulk systems 
can appear. The other one is to introduce a novel bulk 
physical quantity which can be directly related to the 
winding number. This quantity is referred to as chiral 
charge polarization. We show that for 3D class AIII TIs, 
the chiral charge polarization is induced by an applied 
magnetic, which is in analogy with topological magne- 
toelectric effect, and furthermore, the winding number 
appears in its response function (see Eq. (fT6]) below). 

Bulk winding number and magnetoelectric polarization 
in chiral symmetric TIs — We, first, consider the ap- 
proach based on heterostructure systems. To explain 
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our approach in a concrete way, we consider 3D chiral- 
symmetric (CS) TIs, i.e. a class AIII systems. The 
following argument is straightforwardly extended to the 
case of class Dili TSCs. A key idea is to consider a het- 
erostructure system which consists of the 3D CS TI and a 
chiral-symmetry-broken (CSB) trivial insulator with the 
Hamiltonian, as depicted in FIG. 1(a). Here, the triv- 



ial insulator means that P 3 = in the bulk [17j. For 
instance, we can consider the CSB trivial insulator with 
inversion symmetry in the bulk which ensures P 3 = 0. 
To deal with spatially varying hctcrostructurc systems, 
we utilize an adiabatic approach. That is, as long as 
there is a finite energy gap which separates the ground 
state and the first excited states, the interface struc- 
ture can be smoothly deformed to the slowly varying 
one. In the slowly varying structure, the position op- 
erator z in the Hamiltonian can be treated as a parame- 
ter (adiabatic parameter) independent of 3D momentum 
fc, which parametrizes the spatial inhomogeneity of the 
heterostructure. Then, the magnetoelectric polarization 
P 3 {z) is constructed from the adiabatic Hamiltonian of 
the heterostructure system, H (fc, z). -ff (fc, z) interpolates 
between the bulk Hamiltonian of the CS TI, H(k), and 
that of the CSB trivial insulator, HcsB(k), when z is 
varied; i.e. H(k, z) = H(k) when z is a point in the 
bulk of the CS TI, and H(k, z) = H CSB (k) when z in the 
bulk of the CSB trivial insulator. The adiabatic approach 
was exploited before to derive electromagnetic responses 
of the TRI Z 2 TIs from the axion field theory jH [Hj]. 
Our strategy is to extend the adiabatic argument for the 
Z2 non-triviality to the Z nontrivial electromagnetic re- 
sponses. We, first consider the quantum anomalous Hall 
effect. We note that in the heterostructure junction sys- 
tem, the anomalous Hall effect caused by surface Dirac 
fermions is obtained by integrating z-direction under a 
^-independent electromagnetic field, 
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Here zq (zi) is a point in the CS TI (CSB trivial in- 
sulator), and C is a path of ^-integral. Hence the Hall 
conductivity is given by, 
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There are two important remarks. First, although the 
magnetoelectric polarization P 3 is gauge-invariant only 
for mod 1, the line integral of a small difference oiP 3 {z) is 
fully gauge-invariant. Second, cth is determined not only 
by the bulk magnetoelectric polarization at the point zq 
and that at the point z\ , but also by a homologous equiv- 
alence class of the path C. This means that cth depends 
on the microscopic structure of the interface, and is not 
protected solely by the bulk topology. The concrete path 



C is determined by the signs of the mass gaps of Dirac 
fermions on the surface of the CS TI. In our system, the 
mass gaps are generated by the chiral-symmetry breaking 
field induced by the CSB trivial insulator at the surface 
[lit . Here, we consider the case that the sign of the chiral- 
symmetry breaking field, and hence, that of the induced 
mass gaps are uniform on the interface between the CS TI 
and the CSB trivial insulator. More precisely, the Hamil- 
tonian -ffcsB(fc) satisfying this condition is generally ex- 
pressed as ifcsB(fc) = Ho{k) + a(fc)r where a(k) > 
(or < 0) for any fc, and T is the chiral symmetry opera- 
tor mentioned before, and H (k) does not generate mass 
gaps of the surface Dirac fermions. It is noted that T 
itself plays the role of a chiral-symmetry breaking field. 
Then the winding number A3 for H (fc) and the magneto- 
electric polarization P 3 (z) for the adiabatic Hamiltonian 
H{k, z) satisfies the following relation, 
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This relation is one of our central new findings. We 
present a sketch of the proof of Eq.© below. The de- 
tails are given in the supplemental materials (20j . Since 
the magnetoelectric polarization of the CS TI with the 
winding number A3 is P 3 = ^ mod 1, the value of 
P 3 (z = zq) is fixed by chiral symmetry. Also, as men- 
tioned above, the value of P 3 (z — z\) in the CSB trivial 
insulator is fixed to be zero. Due to these fixed bound- 
ary values, J z 1 c dP 3 (z) is adiabatically invariant, which 
means that this quantity is not changed unless the en- 
ergy gap between the ground state and first excited state 
closes. Thus, we can deform H to a flat band sys- 
tem : H 2 = 1. Note that iJcsB(fc) mentioned above 
is deformed to -ffcsB = ±r without closing the energy 
gap. Then, the adiabatically equivalent class of flat band 
Hamiltonian H is given by 



H(k,0) = cos 6Q(k) ± sin 9T, 



(0) 



where Q(k) = 1 — 2P(k), and P(k) is a projection to the 
occupied bands of H(k), and monotonically changes 
from # = 0to# = 7r/2, asz changes from z to z%. It 
is straightforward to show f* 1 c dP 3 (z) = J Q 2 dP 3 {9) = 
i^ 2 - [2^], where P 3 (9) is the magnetoelectric polarization 
ofH(k,9) defined by ©. 

Using Eqs.© and (@| together, we can readily obtain 
the remarkable result that the winding number A3 can 
appear in the quantized Hall conductivity for the het- 
erostructure system depicted in FIG. 1(a), 
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Hence, the Z non-triviality of CS TIs can be detected 
experimentally in this electromagnetic response. 
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FIG. 1: Heterostructure composed of a CS TI (or a CS TSC) 
and trivial insulators (or superconductors) with Hamiltonian 
H = ±r. In (b) and (c), we assume CS TI (TSC) is coated 
by a CSB trivial insulator (SC) so that the interface structure 
has a finite energy gap anywhere. 



We can also apply the formula (JSJ) to the investigation 
on topological magnctoclectric effects which are char- 
acterized not by the Z2 invariant, but by the Z invari- 
ant N 3 . Let us consider the heterostructure system de- 
picted in FIG. 1(b) and (c), which consists of a cylin- 
drical CS TI with its surface coated by a CSB trivial 
insulators. From Eq. ([5]), the magnetoelectric polariza- 
tion of the CS TI coated by the CSB trivial insulator 
is given by P 3 = P 3 (z = z{) - Q c dP 3 {z) = T^, 



which leads the magnetoelectric effect, P 
andM = -^P 3 £, i i-e., 
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Here, E and B are an electric field and a magnetic field 
applied parallel to the axis of the cylinder. The wind- 
ing number successfully appears in the above magne- 
toelectric responses. It is noted that if the system is 
extended without open boundaries and possesses trans- 
lational symmetry, magnetoelectric polarization ([5J is 
gauge-dependent under large gauge transformation so 
that P 3 n- P 3 + n where n is integer. P 3 = T^f m 
Eqs. © and © implies that the particular choice of 
the configuration of the heterostructure as depicted in 
FIG. 1(b) and (c) corresponds to the particular choice of 
the gauge that can extract the winding number of the TI 
in the heterostructure system. 

Case of TRI TSCs — The above argument is also ap- 
plicable to class Dili TRI TSCs in three dimensions. In 
the case of TSCs with spin-triplet pairing, since both 
charge and spin are not conserved, it is difficult to detect 
the topological character in electromagnetic responses. 
However, instead, thermal responses can be a good probe 
for the topological nontriviality, because surface Majo- 
rana fermions still preserve energy. An effective low en- 
ergy theory for the thermal responses of TSCs is the 
gravitational axion field theory described by the action 
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where E g is a gravitoelectric field which play the same 
role as temperature gradient -VT/T, and B g is a gravit- 
omagnetic field, which is in analogy with a magnetic field 
of electromagnetism, and v is the fermi velocity. Because 
of time- reversal symmetry, P3 in the above action (|10j) . 
takes only two values, i.e. or 1/2, implying the Z2 non- 
triviality, and hence Eq. (fTCTj) is an incomplete description 
for the Z nontrivial TSCs. However, as in the case of 
class AIII TIs discussed above, the winding number N 3 
can be detected as thermal responses in a certain class of 
heterostructure system. In the case of class AIII TIs, an 
important role is played by the chiral-symmetry-breaking 
field T. Similarly, also in the case of TRI TSCs, the 
winding number appears in the heterostructure system 
composed of a TRI TSC and a trivial phase with broken 
chiral-symmetry For a p-wave TSC which is realized 
in 3 He, Cua:Bi2Se3 and Li2Pt3B, the chiral-symmetry- 
breaking field is nothing but an s-wave pairing gap with 
broken time reversal symmetry. This is easily seen from 
the fact that the s-wave pairing term of the Hamilto- 
nian is expressed as HeA s T y o- y + ImA s T where A s is the 
s-wave gap, and r M (a^) is the Pauli matrix for particle- 
hole (spin) space, and the chiral symmetry operator T is 
expressed as T = r x a y . When the imaginary part of A s is 
nonzero, this term breaks chiral-symmetry Thus, Eq.Q 
is applicable for the heterostructure system composed of 
a TRI TSC and a trivial s-wave SC with broken time- 
reversal symmetry, as long as the real part of the s-wave 
gap does not yield gap-closing. For the system depicted 
in FIGs.l(a), (b) and (c) the quantum anomalous ther- 
mal Hall effect and the topological gravitomagctoelectric 
effects associated with the winding number are realized. 
Combining the gravitational axion field theory (fTU)) and 
the relation (J5]) , we obtain the quantum anomalous ther- 
mal Hall conductivity, 
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realized for the system shown in FIG. 1(a). This result 
essentially coincides with that obtained by Wang ct al. 
from the argument based on surface Majorana fermions 
fl(| . We can also obtain the gravitomagnetoelectric ef- 
fects. 
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realized for the system shown in FIG. 1(b) and (c). 
Eq. (|12[) implies that circulating energy current flows sur- 
rounding the axis of the cylinder induces the energy (or 
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thermal) polarization, resulting in nonzero temperature 
gradient along the axis. The winding number explicitly 
appears in this thermal response. 

Chiral polarization and the winding number — Hith- 
erto, we have explored the Z topological responses in 
heterostructure junction systems in which the winding 
number successfully emerges as the quantum Hall effect 
and the topological magnetoelectric effect. However, it is 
still desirable to establish a direct connection between the 
winding number and the bulk physical quantities, as in 
the case of the quantum Hall effect in a two-dimensional 
electron gas and Z 2 TIs. We pursue this possibility here. 
For this purpose, we introduce the chiral polarization de- 
fined by, 

P 5 = y c E KI* 5 K), (14) 
' neoccupicd 

where \w n ) is the Wannier function, V c is the unit cell 
volume, and X 5 is the projected chiral position opera- 
tor defined by = PTf^P with f a position operator 
and P the projection to the occupies states. Generally, to 
construct the Wannier function localized exponentially in 
real space, we need the absence of gauge obstruction of 
the Bloch wave function, i.e., vanishing of Chcrn number 
C ij /{2m) = j BZ d 3 k/(2Tr) 3 trJ 7 ij = @|- In chiral sym- 
metric systems, the Chern numbers are zeroes [29(, 
and hence the exponentially localized Wannier functions 
are always well defined. Eq. (|14[) is similar to charge 
polarization, but an important difference is that the chi- 
ral symmetry operator T is inserted in (|14[) . For the class 
AIII TIs and the class BDI TIs with two sub-lattice struc- 
tures, .P 5 represents a difference of charge polarization 
between two sub-lattices. It is noted that in contrast to 
charge polarization which depends on the choice of gauge, 
P 5 is gauge- invariant, since the gauge ambiguity cancels 
out between the two sub-lattice contributions [2(|. As 
will be shown below, _P 5 is a key bulk quantity which 
can be related to the winding number. Actually, in the 
case of one-dimensional (ID) systems, P 5 is expressed by 
the ID winding number N\ as [2oT |. 



For instance, for the ID BDI class TIs such as the Su- 
Schricffcr-Hccgcr model of polyacetylene, Eq. fTS")) repre- 
sents fractional charges which appear at open edges of 
the system. Eq. (|15p is derived from non-trivial algebraic 
properties satisfied by X 5 , which can be regarded as a 
generalization of the commutation relation of the pro- 
jected position operator PrP [20J. In the 3D case, this 
algebra also yields an interesting result that the winding 
number A3 is expressed by the Nambu three bracket of 
X 5 [i(| [5o|], which recently attracts much attentions in 
connection with the density algebra in 3D TIs [3l|, [32| ■ 
However, we have not yet succeeded to relate the Nambu 



bracket to any physical quantities in condensed matter 
systems. Thus, we here take a different approach for 
the 3D case. In fact, in the case of 3D AIII TIs, on 
the assumption that the occupied and unoccupied Wan- 
nier states satisfy the chiral symmetry \wn) = T \ w n ) 
(n G unoccupied, n € occupied ), a more remarkable and 
useful relation between P 5 and the winding number A3 
can be derived; P 5 can be induced by an applied mag- 
netic field, in analogy with the topological magnetoelec- 
tric effect, and furthermore, A3 appears in the response 
function. From the first-order perturbative calculation 
with respect to a magnetic field, we obtain [20l ]. 

P 5 = -£- c N S B. (16) 

Thus, the winding number can be detected as the chiral 
polarization induced by a magnetic field. This is another 
main result of this paper. 

It is expected that an analogous effect may be real- 
ized in 3D TRI TSCs. In the case of TSCs, to explore 
topological characters, we need to consider thermal re- 
sponses, instead of electromagnetic ones. However, we 
have not yet succeeded to obtain thermal analogue of 
Eq. tfTB")) . Furthermore, it is highly non-trivial what P 5 
means for the case of superconductors. These are impor- 
tant open issues which should be addressed in the near 
future. 

Conclusion - We have clarified that the Z non- 
triviality of 3D TRI TSCs and TIs with sub-lattice sym- 
metry can appear in electromagnetic and thermal re- 
sponses of heterostructure systems which consist of the 
TSCs or TIs and CSB trivial s-wave superconductors or 
band insulators. We have also established the relation 
between the bulk winding number and the bulk chiral 
polarization, which may be utilized for experimental de- 
tection of the Z non-triviality. 
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SUPPLEMENTAL MATERIAL 
The derivation of Eq. (5) 

In this section, we present a detailed derivation of Eq. (5). We calculate the continuous change of the magnetoelectric 
polarization Ps(z) between the chiral-symmetric topological insulator H(k) with winding number N3 and chiral- 
symmetry-broken trivial insulator HcsB(k) = Ho(k) + a(k)T with a trivial magnetoelectric polarization P3 = 0. As 
noted in the main text, we assume that a(k) > (or < 0) for any k, and H {k) does not generate mass gaps. Thus 
all of the signs of the mass gaps of the surface Dirac fermions are determined by the sign of a(k). Since the value 
of f** c dPz{z) is adiabatically protected against smooth deformation of surface structure which does not close the 
energy gap, it is sufficient for our purpose to consider a fiat band system, and assume, without loss of generality, that 
the spatial inhomogeneity of the heterostructure system is sufficiently slow, allowing the semiclassical treatment of 
the spatially varying parameter. Then, the Hamiltonian of our system is expressed in the form of Eq. (6) in the main 
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text, 



H(k,9) = cos 9Q(k) ± sin 6T, 



(SI) 



where Q{k) is the "Q-function" defined by Q(k) = 1 — 2P(k) with P(k) the projection to occupied bands of H{k). 
The continuous change of the magnetoelectric polarization P 3 {z) is given by, 



dP 3 (z) 



z ,C 



dP 3 (9) 

de 



8tt 2 



trJ" 2 (fc,6»), 



(S2) 



=0 JBZ 



where P 3 (9) is the magnetoelectric polarization of H(k, 9), T = dA + A A A is the Berry curvature, and A(k, 9) = 
i (u(k, 9) | du(k, 9)) is the Berry connection for the occupied states \u(k, 9)} of the semiclassical Hamiltonian H(k, 9). 

Note that the Q-function of H (k, 9) is equivalent to H (k, 9) itself as shown below. The occupied states \u(k, 9)) of 
H(k, 9) are given by 



\u(k, 9)) = cos - |«(fc)} T sin -T \u(k)} 



(S3) 



with \u(k)) the occupied states of Q(k), Q(k) \u(k)) 
Q(k,9) = 1 - 2P(k,9) 

= 1-2 ]T (\u(k,e))(u(k,6)\) 



\u(k)). Then the Q-function of H(k, 9) is 



unoccupied 

cos9 ( T \ u 

unoccupied 

cos(9Q(fc)±sin6T 



(fc)) (u(k)\T - \u(k)) (u(fc)|) ± sin 6> ^ rrKfc))(u(fc)| + |u(fc))(u(fe)|r 



unoccupied 



(S4) 



Here, we have used the relation TP(A;) + P(fc)T = T obtained from the chiral symmetry. 

We now calculated the right-hand side of (|S2[) . Generally, the Chern form trJ 7 ™ can be written in terms of the 
gauge invariant projection In the case of n = 2, 



tr.F 2 (fc,0) = -tr 



P(fc,0) (dP(fc,l 



or equivalently, with the use of <iPP + PdP = dP and eiQ = — 2<iP, we have, 



1 



trJ^(M) = 55 tr 



Q(k,0) (dQ(k, 



(S5) 



(S6) 



Dividing the external differential d into (d k ,de) = (dk x , dfcj/, d/t. , de), and using Eq. (|S4[) . we rewrite the right-hand 
side of Eq.([S6} as, 



tr 



Q(dQ 



= tr 
= 4 tr 
= 4 tr 
= ±4 tr 



Qj(dfcQ) 3 Ad g Q+ (dkO) A d e Q A d k Q + d k Q A d g Q A (d k Q^j + d e Q A (d k Q 
Q(dhQ) AdeQ 

( cos 9Q{k) ± sin6T) (cos6W fc Q(fc)) 3 A ( - sin 0Q(fc) ± cosdv) d6 



~\ 2 



~\ 2 



_\ 3 



(S7) 



rQ(fe)(d fc Q(fe)) i 



A cos 3 9d9. 



Here, at the second line of cq. (|S7[) . we used a cyclicity of the trace and d k 

used rQ(fc) + Q(k)T = followed from the chiral symmetry. For the basis in which T is represented as T = 



Qd k Q, and at the forth line, we 
1 



Q(k) is expressed as Q(k) = 



q(k) 
gf(fe) 



8tt 2 



trJ" 2 (fc ; 



-lj' 

with a unitary matrix q(k). Then, from Eqs. (|S6j) and (|S7|I . we have 
1 



"32tt 



rtr y (k)d k q(k)] A cos 3 0d0. 



(S8) 
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Hence the change of the magnetoelectric polarization (|S2[) is given by 

P dP 3 (z) = / tr [ q Hk)d kq (k)} 3 f 

Jzn.C 6Z7T JBZ J0 = 



cos 3 



= ±^-L J Bz tr[ q Hk)d k q(k)] 3 (S9) 

where -/V3 = Jgz^ T [rf (tydkOlk)] is the winding number characterizing the ground state topology of the chiral 
symmetric topological insulator [2j. 



The winding number in arbitrary odd spacial dimensions 

The winding number A^n+i of the chiral symmetric topological insulator (superconductor) in 2n + 1 spacial di- 
mensions characterize the homotopy of the map from BZ 9 fc n- q(k) in the unitary group U(m), where q(k) is the 

off-diagonal part of Q(k) on the basis such that T is represented as T = ^ . N2 n +i is given by 

^n+l = ~ , . n+10 L , ,s„ / tr [q\k)dq(k)] 2n+1 . (S10) 



(27ri)" +1 2"(2n+ 1)!! 



Eq. (|S10[) is rewritten as 



^2n+l = 7 r-n TT7 777 / tr T \Q(k)dQ(k) + 

(Sll) 

tr rg(fc) [dg(fc)] 2n+1 . 



(-i) n 



(27ri) n+1 2 n+1 (2n+ 1)!! 
Here, we used Q(k)dQ(k)Q(k) = -dQ(fc) from Q 2 (fc) = 1. 



The derivation of Eq. (15) and some algebraic properties of projected chiral position operator 



In this section, we will derive Eq. (15). For this purpose, we, first, explain some important features of the 
projected chiral position operator X^ which hold both in ID and 3D systems. X^ has some similarity with the 



projected position operator X^ = Pf u P where P = Z^neoccupiedSfeeBZ l^h) (4>nk\ is a projector on the occupied 
bands. It is useful for the following argument to summarize some basic properties of X^ here. On the basis of Bloch 



states \4> n k), is represented as {4> nk 
with \rik) 
Xfj, yields 



X„ 



iS 



a 

nm dk u 



A nm ^(k)) 8k,k', where A nm ,n{k) 



»k 



(>ni k 

dk u 



-ik-r 



ifc) is the Berry connection. Then the non-commutativity of the projected position operator 



X^, X v 



k> J = iF nm ^ v {k)8 k ^ , where J>„(fc) = d u A v - d v A u + [A^,A U ] is the non-Abelian 
Berry curvature arising in multi-band systems. Under a gauge transformation of the occupied Bloch states Vl'(fe) = 
{\n k ) , \rrik) , ■ ■ ■ } n m ... gocc l— > ^(k)U(k) with a unitary matrix U(k), the Berry curvature T is transformed as T 1— > 
U'TU, hence a nontrivial gauge dependence do not exist. 

We now apply a similar argument to the projected chiral position operator defined by 



X* := PTr^P. 



(S12) 



The representation of on the basis of Bloch states is given by 



'nk 



x; 



= { n k 



—ih 



n k 



dk,. 



d 



->k,k' 



Ah 



dk' 



m kl 



(S13) 



8 



Note that {4> n k 1 1 
Wc denote X^(k) := 



vanishes when both \<j> n k) and \4> m k') are occupied states, because of the chiral symmetry. 

or Xf t (k) := (k)Td tl 'i'(k) in matrix representation, Hereafter, X^(k) is 



i(n k 



dk u 



referred to as the projected chiral position. Under a gauge transformation ^(k) n- ty(k)U(k), Xf t (k) transforms as 



Xl(k) i y it/'t(fc)\Brt(fe)r5 M {^(k)U(k)} 

= iU f (k)^(k)Td^(k)U(k) + iC/ t (fc)* t (fc)r*(fc)a AI L/(fc) 
= U\k)Xl{k)U{k). 



(S14) 



Thus, Here, we used (rife | T | mj~) = for the n, m occupied bands. Hence, in the chiral symmetric systems, the 
projected chiral position X^{k) is gauge invariant in the same way as the Berry curvature ^^(fc), in sharp contrast 
to the projected position which is gauge- dependent. Furthermore, the commutator of the projected chiral position 
yields the Berry curvature, 



-^Td^^Td,^ - (p O v) 
(1 - ** f ) d v ^ - (fi o v) 



(S15) 



Here wc used 



E i n ) ( n i 

n^occpied 



E 



\n) (n\ 



nGunoccupied 



£ [ 



n) (n| +T|n) (n|T 



unoccupied 



+ r** t r. 



(S16) 



We also remark that the winding number is written as the integral of the (2n + 1) -bracket of the projected chiral 
position over the Brillouin zone : 



N- 



2n+l 



7r n + 1 (2n+ 1)!! JBZ 
-i n e„ 



tr 



A l> A 2i' 



y-5 



C MlA t 2---M2n + l 

Tr n+1 {2n + 1)!! 



j2n+l 



fc tr 



BZ 



Ai2n + 1 



(S17) 



where [Xl, X 2 , • • • , X 2n +i] = ^i^-^n+iX^X^ ■ ■ ■ X^ 2n+1 is the (2n + lV-bracket, and /Zj run over 1, 2, • ■ • , 2n + 1. 
Note that when n = 1 (i.e. A^), [Xi, X2, X3] is the Nambu bracket [3j. This expression of the winding num- 
ber is analogous to the relation between the Chern number and non-commutativity of the projected position 

operator. Eq. (|S17|) directly follows from (|Sllj) . The matrix element of dQ = d (l — 2 Ylneoccupied l n ) ( n l) = 

~^Y1 unoccupied ( M n ) ( n \ + \ n ) (^ n l) 011 the basis of the occupied states \n) and unoccupied states |n) = T |n) are 



(n | dQ | n') = 0, 
(n\dQ\ n) = 0, 
{n\dQ\n') =2{n\V\ dn') , 
(n | dQ | n') = — 2 (n | T | dn') . 

Hence dQ is written as the matrix form in the space spanned by the occupied and unoccupied states. 



(S18) 



dQ 



(n\dQ\ n') (n\dQ\ n')\ f 2 (n | T | dn') 

(n I dQ I n'} (n | dQ \ n')J ~ \-2 (n [ T | dn') y ( S19 ) 
= -2ir 2 (n | T | dn') = -2r 2 l r 5 m ^ dfc p 

where r = (ri, r 2 , T3) is the Pauli matrices in the occupied-unoccupied space. In the same way, T and Q is written as 

5 n 



[n\T\n') (n\T\n>) 
n\T\n>) (n\T\n>) 







Tl5 n 



(S20) 



\n\Q\n') (n\Q\n') 
[n\Q\n') (n\Q\n') 



-T 3 Sr, 



(S21) 
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By using (|S19|) . (|S20[) and (|S21[) . we can write the winding number (|S as 



2n+l 



(-iy 



(2 7 rz)"+ 1 2"+ 1 (2n+ 1)!! j 

(-1)" f 
(2 7 ri)"+ 1 2™+ 1 (2n+ 1)!! J BZ 

(«)«(2n+l)!!7 BZ 



tr rQ(fe) [dQ(fc)] 



2n+l 



tr 



n(-T 3 ) (-275^ d* 



tr 



A M1 A M2 ••• A M2„ + l 



C MlM2---M2n + l 



(7ri) n (2n + 1)!! 



d 2 " +1 fe tr 



A /*1 A M2 ' ' ' A M2„ + 1 



Thus, we arrive at (|S17|) . 

In the case of the one dimensional systems, the winding number 



iVi 



1 



tr X 5 (k) 



(S22) 



(S23) 



T j BZ 

is directly related to the chiral polarization defined by (14). For the Wannier function localized at the site R, (14) is 



P' 



V C N C 



E E 

unoccupied k,k f 



e lkR Yxe- lWR 



V C N C 



e 

27 



unoccupied k,k' 

dk tr X 5 (k) 

BZ 



(S24) 



TVie 
~ 2 ' 

where iV c is the number of unit cell. Thus we obtain Eq.(15). 

To close this section, we would like to comment on an implication of (|S17|) to the 3D case. In this case, from 
Eq. (|S17[) . the winding number N3 is expressed in terms of the Nambu bracket. This implies that N3 may be related 
to a physical quantity described by the Nambu dynamics 0]. However, we do not know any examples in condensed 
matter systems which are described by the Nambu mechanics, and also quantum version of the Nambu mechanics is 
not well understood. We have not yet succeeded to obtain any insight from the study in this direction. Thus, we 
consider a different approach to relate the 3D winding number N3 to a physical quantity, as explained in the main 
text and the following section. 



The derivation of Eq. (16) 

In this section, we derive Eq.(16) on the assumption that the occupied and unoccupied Wannier states satisfy the 
chiral symmetry \iUnR) — T\w n R) (n <= unoccupied, n € occupied ). This is equivalent to the gauge fixing condition 
\4>nk) = r \4>nk) for unoccupied and occupied Bloch states. Our derivation of Eq.(16) is based on the perturbation 
formalism developed by Kita-Arai for the Wannier function under the magnetic field. Generally, to construct the 
exponentially localized Wannier function, we need the absence of gauge obstruction of the Bloch wave function, i.e., 
vanishing of Chern number Cij/(2ni) = J BZ d 3 fc/(27r) 3 trJ r y = 0. In chiral symmetric systems, the Chern numbers 
Cij are zeroes, and hence the exponentially localized Wannier functions are always well defined. 

The chiral charge polarization in the case with no magnetic field introduced in the main text is 

P (B=o) = y E ( WnR I rf I WnR ) ' (S25) 

c n£occ 

where \w n n) is the Wannier function localized at a site R constructed from the n-th occupied band, 

\w nR ) = — j= e ~ lk R \M . (S26) 
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where N c is the number of unit cells. Now, we will calculate the first order perturbative corrections to P 5 with respect 
to an applied uniform magnetic field B. First, we introduce modified Wannicr states \w' nR ) defined by [f| 

<r(t) = e u ^w nR (r), (S27) 

where I rR is the Peierls phase 

I rR = Y c j R dr ' ■ A (r') (S28) 

with dr' the straight line path from R to r. The modified Wannier states \w' nR ) form a complete set, though they 
are not orthonormal for the case with a finite magnetic field. To orthonormalize them, we need to include corrections 
from other sites and other bands to \w' nR ). Then, the orthonormal modified Wannier function is expressed as, 

\<PnR) = E \ W 'n'R') Sn'R',nR, (S29) 
n'R 

for which (<p n R\ip n ' R 1 ) = 5 nn '8 RR i is satisfied. Thus, the chiral charge polarization under an applied magnetic field is 
given by 

p5 = y E & nR I rf I ■ (S30) 
From Ref. in the cases of a uniform magnetic field, Su'R^nR * s expanded up to the first order in B: 

il i 

Sn'R.nR = 5 nn ,S RR , - ^-B^f—— £ e *(*'-«) (d^dwk) , (S31) 

C c k 

where \nk) = e~ lk ' r \4> n k)- We denote \<p n a) = \ w ' nR ) + \^ w nR)- Then, the correction term of the chiral polarization 
is 

SP 5 = y E ^ w 'nR I r * I <r) - Kfl I Tr I w nR )} + Y WnR I r * I Sw 'nR.) + c - c -l ■ (S32) 

c jiGocc c n£Locc 

The first term vanishes since the Peierls phases of \w' nR ) are canceled out in (w' nR | Tf \ w' nR ). The second term is 
recast into the following form with the use of (|S31[) . 

«r I Tr I Sw' nR ) = E Mn I ^ I <'R') ^ E e* (R '- R) (S,n' fc |ftn*) 

n'R' c k (S33) 

n'R' k 

Here, we omitted the Peierls phases since they are higher order corrections. On the other hand, the factor 
(wnR | Tf | Wn'R') is expressed as, 



(w nR | Tf | w n , R ,) = -L e lk R ((jink | Tf \ cp n , k ') e~ ik : ' R ' 

c k,k' 

= W E eik R {" f ^7 (WK> <>*.*') + 1 W\Vn' k ) S k ,k' } e- lk ' R ' (S34) 
= RSnn'S R ,R + yJ2 e lk < R - R 'h (nk\T\Vn' k ) , 



N c 

k 



where n is a label for an unoccupied state and we fixed the gauge of \4> n ) satisfying \<fi n ) = V \cf) n ) with \<$> n ) an occupied 
state. From (|S33[) and (|S34j) , we have, 



«r I Tr | <5< H ) = -^LB&fl— E 



fc 



(S35) 
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The first term vanishes because it is the total derivative : e^i (djrik\Y\dink) = eijidj {(nk\Y\dirik)} ■ Hence we obtain, 



.-,2 



5P 5 = JL-B^j^y E M r l V <> %n' k \d ink ) + c.c. 

c neoccn'k /g 36 N 

r d 3 k 



f d 3 k 

e yj i J -— gtr [^Yd^d^d^ + ^d^d^Yd^] 



2hc 

where VP = {\ni) , 1 77,2) , • • • } n .c OC c an ^ we om it the parameter k. 

From this expression, we find that ay = j^- is proportional to Sij. For instance, we consider an off-diagonal 
component, 

d 3 k 

y ~Wc 

t Yk (S37) 

= / - — -tr \&Td x Vd z &d x V - &rd x Vd x &d g V + &d x Vd z &Td x V - &d x Vd x &Td x V] . 

2nc J (2ir) 3 1 J 

We see that on the interior of the trace in Eq. (|S37l) . 

(1st) + (4th) = tr [&Yd x yd z &d x y - ¥ d x ^d x ^Yd z ^] 

= tr [-d x Vd x &Y (*0,tft + d M V&)] (S38) 
= tr [-d x *d x &Yd z P] , 



and 



(2nd) + (3rd) = tr [-^Yd x ^d x ^d z ^ + ¥d x ^d z ^Yd x ^f] 

= tr [-d x m x ¥ (d z ^ + T] (S39) 

= tr [-5 B *flS B *t fl# pr] . 



Thus, we can show 



(S40) 



where, P = X^neocc l n ) ( n l * s a projector to occupied bands, and we have used dW^ + Vd^ = dP, and Yd z P+d z PY = 
d z Y = 0. Hence, 



e 2 /" 

487r3fie / tr [* t rd*d* t d1' + * t d1'd* t rd*] 
^ / tr [* +rMtd *] ' 



(S41) 



The integrand of this expression is gauge-invariant under gauge transformation \& — > i&U with a unitary matrix U 
except for total derivative. This is seen as follows. 

tr [tftlWdtftdtf] 
i-> tr [J7 t * t rrf(*C/)d([/ t ^ t )d(\E'i7)] 

= tr [&YdV (d¥d^ + UdU^dtif + VUVUdtf - UdU^UdU^)] 

= tr [¥Yd^d¥d^] + tr [¥ d^^YdWdU^ + ^Yd^^dWdU^ + ^YdWdUdrf] (S42) 
= tr [^Yd^d^d^] + tr [-d^ (** t r + T** 1 ) d^UdU^ + &YdVdUdU f ] 
= tr [^Yd^d^dV] + tr [-dtftlWUeff/t + ^Yd^fdUdU^] 
= tr [* t rtZv&d* 1 'd*] - d tr [ftrdWeff/ 1 ] , 

where, we have used V&Y + r#*t = £ n£occ (N ("I r + r | n) (n\) = Y. Hence the integral over the Brillouin zone is 
gauge- invariant. Now, we choose the basis, 
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where, q is the off diagonal part of Q-function Q 



( 



q 




;) 



. Then 



a 




I 



tT(q^dq) 3 



(S44) 



where iV 3 = -^h^ I ^iq^dq) 3 is the winding number. Thus we arrive at Eq. (16). 

Note that the magnetic induced chiral polarization represented by the first line of (|S36[) generally depends on the 
gauge of the unoccupied Bloch states \4>nk), so the gauge fixing condition \4>nh) = T \4> n k) is crucial for our results. 
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